Classical theory of second-harmonic generation from magnetic metamaterials 
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Strong second-harmonic generation has recently been experimentally observed from metamateri- 
als consisting of periodic arrays of metal split ring resonators with an effective negative magnetic 
permeability [Science, 313, 502 (2006)]. Because of the appearance of the fundamental plasmonic 
resonance, both electric and magnetic fields are strongly localized inside gap of the resonator, which 
in turn greatly enhances the conversion efficiency of second-harmonic generation. To explore the 
underlying physical mechanism, a classical model derived from microscopic theory is employed here. 
The quasi-free electrons inside the metal are approximated as a classical Coulomb-interacting elec- 
tron gas, and their motion under the excitation of an external electromagnetic field is described by 
the cold-plasma wave equations. Through numerical simulations, it is demonstrated that the micro- 
scopic theory includes the dominant physical mechanisms both qualitatively and quantitatively, ft is 
shown that, in contrast to previous interpretations, the dominant source term for second-harmonic 
generation results from the convective derivative of the continuous electron current. 

PACS numbers: 42.70.-a, 52.35.Mw 



Strong second-harmonic generation (SHG) from mag- 
netic matematerials (MMs) has recently been experimen- 
tally observed [l|, 0, Q- ^^^ to the periodic arrange- 
ment of artificial metallic subwavelength structures, such 
MMs possess plasmonic resonances as well as frequency 
regions with effective negative magnetic permeabilities 
0) S M, 0]- III particular, under the excitation with 
an external electromagnetic (EM) wave, quasi-free elec- 
trons inside MMs are coherently driven by the Lorentz 
force and form a circulating current. The current in 
turn generates an induced magnetic field pointing op- 
posite to the original one. In a first explanation of the 
observed intensity dependence of the SHG on sample ge- 
ometry and polarization direction of the exciting field, 
it has been assumed that the magnetic component of 
the Lorentz force provides the dominative second-order 
source [l|, 01 • However, in a theoretical description of a 
classical Coulomb-interacting electron fluid, three differ- 
ent nonlinear source terms exist which can potentially 
contribute to the observed second-harmonic (SH) sig- 
nal. In this Letter, we derive and numerically solve 
the corresponding equations and study the underlying 
physical mechanism of SHG theoretically. We show that 
our model gives good agreement with the experimentally 
observed SH strengths and thus catches the dominant 
physical mechanism. A switch-off analysis furthermore 
demonstrates that the dominant SH source can be traced 
back to the convective derivative of the continuum gen- 
eralization of Newton's equation. 

The motion of quasi-free electrons inside a metal is de- 
scribed classically. Quantum mechanical Coulomb corre- 
lations and exchange contributions are thus missing while 
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the classical Coulomb interaction (i.e. the Hartree term 
in a quantum mechanical derivation) is fully included [8| . 
Furthermore, the Coulomb scattering due to higher-order 
quantum corrections is phenomenologically described via 
an inverse decay time 7 = 1/r [y,[iOjllli- The electrons 
inside the metal are approximated as a continuous fluid 
and macroscopic quantities such as electron number den- 
sity Ue and electron velocity Ug are therefore continuous 
functions of position r and time t [i2jll3]- The fluid satis- 
fies the conservation of particle density which is described 
by the continuity equation 
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Under the excitation of an external EM wave with elec- 
tric field E and magnetic field B, the dynamics of the 
electron fluid is governed by the continuum generaliza- 
tion of Newton's equation. 
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Here, me and — e are the electron mass and charge, re- 
spectively. The so-called convective or material deriva- 
tive of Ue on the left-hand side is a result of the descrip- 
tion of electrons in terms of a continuous density and 
can be formally derived from the quantum mechanical 
Wigner distribution [a|. More intuitively, it can be un- 
derstood as the time derivative of an electron taken with 
respect to a coordinate system which is itself moving with 
velocity Ue(r(t),i), given by J12l.ll3l| 
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To describe the interaction between the electron fluid 
and the EM field self-consistcntly, the electron dynam- 
ics have to be coupled to Maxwell's equations. Conse- 
quently, it is convenient to express the charge density 



and the current density, 

p{r,t) =e(nj(r)-?ie(r,i)), 

j(i">0 = (-e)ne(r,t)ue(r,t), 



(4) 
(5) 



in terms of the number density and the velocity field. 
Here, the ions are assumed to have an infinite mass such 
that only the electrons contribute to the electric current 
and the ion number density rii is time-independent. By 
choosing rii equal to a constant inside the metal and zero 
outside, the geometry of the metallic MM is fully de- 
scribed by ni{r). 

Using the definitions from Eqs. (|4]) and ([5]) together 
with the differential Eqs. ([1]) and ([2]), we can derive 
equations for p and j. The total set of these equa- 
tions coupled to the EM fields via Maxwell's equations is 
known in plasma physics as the cold-plasma wave equa- 
tions [12l . Il3| . The equations are 
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In practice, the set of Eqs. ([6|l- ([TT|) is redundant to a 
certain degree. For example, Faraday's law of induction, 
Eq. ([5]) , automatically guarantees a vanishing divergence 
of B for all times if B is initially chosen divergence-free. 
In the same way, Eq. © along with Eq. ^ are equivalent 
to the continuity Eq. pT|) [14]. Hence, only the reduced 
set of Eqs. (I71)-([Tni) has to be solved. 

By expanding the above equations in terms of the peak 
electric-field amplitude of the illuminating EM pulse, we 
can obtain two sets of equations [l^ [la, [i3 • The first 
set describes the linear-optical (first-order) response of 
a classical Drude-type metal whose permittivity is given 
by 



e{uj) = 1.0- 
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The second set, taking the linear fields (Ei, Bi and ji) 
as sources, captures the second-order nonlinear-optical 
process of the metal. 
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where the bulk plasma frequency tjp(r) is defined by 
ujp{r) — e^ni(r)/eo?7ie [16[, and k represents x, y and z 
coordinate. The form of this set of equations immediately 
suggests: (i) its homogeneous part is exactly identical to 
the linear set such that the propagation of the gener- 
ated SH wave is modified by the Drude response of the 
metal; and (ii) there are three distinct nonlinear sources 
for SHG in the plasma, representing different physical 
mechanisms. The second and third source are the well- 
known electric- and magnetic-component of the Lorentz 
force, respectively. The first source term is a generalized 
divergence originating from the convective time deriva- 
tive of the electron velocity field Ug mentioned above. 

For the numerical solution of the coupled differential 
Eqs. P^ - P3|) and their linear counterparts, a three- 
dimensional finite-difference time-domain (FDTD) al- 
gorithm is applied. According to Yee's discretization 
scheme, the whole computational space is divided into 
millions of cubic grids [18] . Due to the translational sym- 
metry in the xy plane it is sufficient to simulate a single 
unit cell of the periodic MM. Periodic boundary condi- 
tions are employed in the zy-direction while perfectly- 
matched layers are implemented in z-direction perpen- 
dicular to the metallic layer. The incoming electric field 
propagating along the —z direction is generated by a 
total-field/scattering-field technique [ig| . 

A variety of metallic MMs with different configura- 
tions have been studied numerically, and the details of 
these simulations will be discussed in a second publica- 
tion [l7|. From a detailed quantitative comparison be- 
tween numerical simulations and experimental measure- 
ments, it is concluded that our simulations agree very 
well with the experimental observations [ll, 0, S]- ^'^^ 
instance, we have confirmed the symmetry requirement 
that far-field SHG signals polarized along a direction of 
mirror symmetry are forbidden. Thus, a rectangular ar- 
ray of "U" or "T" structures can only show SHG in one 
polarization direction while SHG is absolutely forbidden 
from arrays of centro-symmetrical structures. Thus, our 
theory of the classical Coulomb-interacting electron gas 
catches the physical mechanism of SHG from metallic 
MMs. 

In the following, we focus on arrays of metallic split- 
ring-resonators (SRRs) similar to those studied in Ref.[ll]. 
Those materials are known to possess negative effective 
permeabilities in certain frequency regions and are gener- 
ally referred to as magnetic MMs. A schematic drawing 
of a single unit cell of the periodic SRR array is shown 
in Fig.l. The complete structure then consists of a rect- 
angular lattice of SRRs periodically arranged in the xy 
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FIG. 1: (color), (a) Schematic drawing of a rectangle lattice 
of gold split-ring resonators, (b) There are three layers (gold, 
ITO and glass) along the z direction. 



plane. The cross section of the structure in Fig. (lb) 
shows that the nietalhc layer is supported by a glass sub- 
strate coated with a thin film of indium-tin-oxide (ITO) 
in the z-direction. Furthermore, the SRRs are made of 
gold, and the thicknesses of the gold layer and the ITO 
film on the semi-infinite glass substrate are 25 nm and 5 
nm, respectively. 

To study the dependence of SHG on the different-order 
plasmonic resonances, two different-size SRRs are consid- 
ered. The small one, named as SI, closely matches the 
small-SRR sample fabricated in Ref.[l|, has lattice con- 
stants of the rectangle lattice ax = ay = 305 nm, the 
width of the metallic ring w = 60 nm, Ix = ly = 220 nm 
and 6=115 nm. The big one, referred to as S2, resembles 
SI graphically and takes geometric parameters as follows: 
ttx — 610 nm, ay — 635 nm, Ix — 455 nm, ly = 470 nm, 
w — 135 nm and b — 95 nm. In addition, the permittiv- 
ity of the ITO layer and glass substrate is 3.8 and 2.25, 
respectively. The bulk plasma frequency of gold is taken 
as LUp = 1.367 X lO^^s"^, and the phenomenological col- 
lision frequency 7 = 6.478 x lO^'^s""'^ [a, 0]. The sizes 
are chosen such that both structures have a resonance 
around a wavelength of 1500 nm. For the purposes of 
physical clarit}^S2 is slightly different from the big-SRR 
sample in Ref.[ll], and geometrically identical to SI. We 
will show below that this disagreement does not infiucncc 
our conclusions qualitatively. 

Our main results are plotted in Fig. 2. Firstly, we study 
the linear-optical properties of these magnetic MMs. 
The transmitted and reflected spectra under illuminating 
light with different polarization are calculated and pre- 
sented. Our numerical spectra agree qualitatively and 
quantitatively with the experimental spectra in Ref.[ll]. 
Our main interest is the transmission null around 1502- 
nni wavelength of SI, and we discuss its physical mech- 
anism in detail. As already pointed out, such a trans- 
mission null is a sign of the fundamental plasmonic reso- 
nance (its mode distribution is shown in Fig. (3a)), which 
is induced by the oscillation of circulating currents in- 
side the metallic rings. This plasmonic resonance bears 
an analogy to an inductor-capacitor circuit resonance. It 
is consequently mainly determined by the inductance of 
the loop and the capacitance of the SRR gap. On the 
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FIG. 2: (color). Summary of the linear-optical spectra 
and normalized second-harmonic intensities for two magnetic 
metamaterials, SI and S2. Both transmission (solid) and 
reflection (dashed) spectra are plotted. The wavelength of 
the illuminating fundamental-frequency wave is flxed at 1502 
nm. The metamaterial SI is excited by (a,a') x- and (b,6') y- 
polarized incidence. The metamaterial S2 is excited by (c,c') 
X- and (d,d') j/-polarized incidence. The polarization of the 
generated second-harmonic waves from both metamaterials 
are only allowed in the y direction. 



other hand, the transmission null of S2 results from the 
second-order plasmonic resonance whose mode is charac- 
terized as a single-node standing wave around the SRR 
metallic ring [3, [l3| , as can be found from its near-field 
distribution plotted in Fig. (3b). 

Next we study SHG from these magnetic MMs illumi- 
nated by a fundamental-frequency (FF) wave with wave- 
length of 1502 nm as well as peak electric- field amplitude 
of 2 X 10^ (V/m) (same as that employed in the experi- 
ments [il, 3). In order to describe the energy conversion 
efficiency in this nonlinear-optical process, we define a 
normalized SH intensity. 



|E2(2c^o)|^ 
|Ei(c.o)P 



(16) 



to measure the strength of the far-field negative z- 
propagating SH wave, where luq is the frequency of the in- 
cident FF wave. Because the infinite-thickness glass sub- 
strate is utilized here, the true SH conversion efficiency 
equals 77 times ^^Cgiass UM- Two important conclusions 
can be observed from Fig. 2: 

(1) Because these magnetic MMs possess x-coordinate 
mirror symmetry, only a j/-polarized SH wave can be ex- 




Min 



Max 



FIG. 3: (color). Magnitude (|-Ba;|) distributions of (a) the 
fundamental plasmonic resonance of SI and (b) the second- 
order plasmonic resonance of S2, at a distance of 5 nm above 
the gold film. Dashed lines mark the positions of the gold 
split-ring-resonators. 



cited by either x- or y-polarized incident fields. Such 
an incidence-polarization dependence was already exper- 
imentally observed [ll, 0, Q. Moreover, a universal selec- 
tion rule can be extracted from all our numerical simu- 
lations |17|, that is, a mirror symmetry of the metallic 
MM in one direction completely prohibits SHG in the 
same direction. 

(2) As shown in Fig. (2a'), the strongest SHG signal 
with a normalized SH intensity of 6.7 x 10~^^ is achieved 
due to the excitation of the fundamental plasmonic res- 
onance of SI by a;-polarized incidence. As a reference, 
the corresponding experimental observation is roughly 
2.0 X 10"" [i,|l|. We now normahze this signal to 100% 
and relate all other signals to it. 

For excitation of the second-order plasmonic resonance 
with x-polarized illumination in Fig. (2c'), a relatively sig- 
nificant SH emission is observed from S2 with normalized 
SH intensity of 18%. (An MM which closely matches 
the big-SRR sample in Ref.[l| provides an SH signal of 
7.8% [17|, quite close to the experimental measurement 
4.3%). On the other hand, only weak SH signals with 
0.2% intensities are found from the off-resonance configu- 
rations, as shown in Fig.(2&',2(i'). Therefore, the presence 
of plasmonic resonances and thus circulating currents 
can greatly enhance SHG from magnetic MMs. More- 
over, different-order plasmonic resonances make different 
contributions. The fundamental resonance contributes 
about fourfold to the second-order resonance. It can be 
interpreted by the localization degree of the linear field. 



because such a compression can enhance SHG even with- 
out perfect phase matching [T^, [l^. As shown in Fig. 3, 
the linear electric fields of both resonances are localized 
inside the SRR gap, and the stronger compression ap- 
pears for the fundamental plasmonic resonance. 

Now we try to distinguish the dominant source among 
these nonlinear terms for SI on resonance. By numeri- 
cally switching off these sources one by one, their respec- 
tive contribution can be roughly estimated. The mag- 
netic Lorentz force, j x B, makes a weak contribution 
as small as 2%. Its electric counterpart provides 29%. 
The most significant contribution, 68%, comes from the 
convective derivation of current. The sum of the last 
two terms, consequently, constructively dominates SHG 
emission from magnetic MMs. It is further found that, 
the contribution of the magnetic Lorentz force is neg- 
ligible as long as the structural fundamental plasmonic 
resonance is excited [l7[. Such an observation is quite 
different from the previous interpretations, in which the 
magnetic contribution was emphasized [l| . 

In conclusions, a classical theory of second-harmonic 
generation from magnetic metamaterials is presented in 
this Letter. The quasi-free electrons inside the metallic 
nanostructure are approximated as a classical Coulomb- 
interacting electron gas, and their motion under the ex- 
citation of an external electromagnetic field is described 
by the cold-plasma wave equations. Through numerical 
simulations, it is demonstrated that the classical theory 
includes the dominant physical mechanisms both qual- 
itatively and quantitatively. It is further shown that: 
(1) A mirror symmetry of metamaterial in one direction 
completely prohibits the generation of far-field second- 
harmonic wave in the same direction; (2) The presences of 
structural plasmonic resonances can significantly enhance 
second-harmonic generation; (3) As long as the lowest- 
order plasmonic mode is excited, the magnetic Lorentz 
force makes negligible contribution to second-harmonic 
generation. 
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